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SECTION  1 


INTRODUCTION 


1.1  PURPOSE 

This  report  describes  some  initial  results  in  the  algorithmic  gener- 
ation of  finite  length  pulse  burst  signals  whose  autocorrelation  functions 
are  approximately  impulse-like  in  nature.  The  codes  will  be  restricted  to 
be  real  with  piecewise  constant  amplitude  that  can  either  be  positive  or 
negative  (i.e.,  bi-phase).  Such  codes  find  various  applications  in  radar 
signal  processing  and  digital  communication  systems. 

1 . 2 BACKGROUND 

In  pulsed  radar  systems  where  it  is  required  to  improve  upon  the  prob- 
ability of  target  detection  for  a given  false  alarm  rate,  it  is  necessary 
to  increase  the  transmitted  pulse's  energy  level.  Ideally  the  transmitted 
pulse  should  have  a very  small  time-duration  and  large  amplitude  in  order 
to  achieve  a high  probability  of  detection  while  maintaining  acceptably  good 
range  resolution  characteristics.  Unfortunately,  since  all  radar  transmit- 
ters have  a peak  power  limitation  of  some  sort,  it  is  not  possible  to  arbi- 
trarily increase  the  transmitted  pulse's  amplitude.  To  raise  the  energy  of 
transmission,  it  is  then  necessary  to  increase  the  pulse's  time-duration. 
This  will  typically  have  a detrimental  effect,  however,  on  the  resultant 
range  resolution  and  ambiguity  function  characteristics. 

There  exists  a class  of  techniques,  generally  refered  to  as  pulse  com- 
pression methods,  that  overcome  the  above  mentioned  inherent  pitfall  which 
large  time-width  pulses  are  susceptable  to.  In  particular,  by  appropriate 
signal  processing  (e.g.,  matched  filtering),  it  is  possible  to  achieve  the 
effect  of  a short  time  duration-large  amplitude  pulse  while  actually  trans- 
mitting a relatively  long  time  duration-moderate  amplitude  pulse.  Two  of 
the  most  widely  studied  of  pulse  compression  methods  are  linear  frequency 
modulating  pulse  compression,  and,  pulse  burst  signals  (e.g.,  see  ref.  [6]). 
In  this  report,  we  will  be  concerned  with  presenting  an  algorithmic  method 
for  synthesizing  pulse  burst  signals  whose  associated  autocorrelation  func- 
tion approximates  the  ideal  spike-like  time-behavior  which  is  the  goal  of 
many  radar  system  designs. 
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SECTION  2 


RADAR  SIGNAL  DESCRIPTION 


2.1  BASIC  SIGNAL  MODEL 

In  radar  target  detection  and  tracking  applications,  the  return  signal 
from  a target  is  often  modeled  as  a simple  transformed  version  of  the  trans- 
mitted signal.  If  the  target  is  taken  to  be  a moving  point  in  space,  then 
the  return  signal  r(t)  is  usually  related  to  the  transmitted  signal  s(t) 
based  on  the  following  two  hypotheses. 

(i)  the  return  signal  will  be  a time  delayed  version  of  the 
transmitted  signal  (i.e.,  r(t)  = s(t-i))  where  t is  propor- 
tional to  the  distance  between  the  target  and  the  transmitting 
radar. 

(ii)  there  will  be  imposed  a frequency  shift  of  on  the  trans- 
mitted signal's  spectrum  which  is  proportional  to  the  target's 
radial  velocity  relative  to  the  radar  system  (i.e.,  R(w)=S (w-w^) ) . 

The  quantity  is  refered  to  as  the  Doppler  shift  and  is 
positive  for  targets  moving  toward  the  radar  system. 

The  combined  effect  of  these  two  factors  will  then  characterize  the  return 
signal  from  a point  target  moving  at  a constant  range  rate.  It  is  tacitly 
assumed  in  this  model  that  the  transmitted  signal  is  narrow  band  and  that 
the  target's  range  rate  is  much  less  than  the  speed  of  light. 

To  demonstrate  the  implications  of  this  signal  model,  let  us  consider 
the  important  case  in  which  the  transmitted  signal  is  a modulated  sinusoid, 
that  is 

s ( t)  = a(t)  cos[uot  + 0(t)]  (2.1-1) 

where  a(t)  and  0(t)  are  the  real  envelope  and  phase  modulating  signals,  re- 
spectively, while  a)Q  is  the  associated  carrier  frequency.  Generally,  these 
two  modulating  signals  vary  slowly  in  time  relative  to  the  unmodulated  car- 
rier signal  cos[u)Qt].  According  to  the  return  signal  model  as  specified  by 
the  above  hypotheses,  it  follows  that  the  associated  return  signal  (in  the 


absence  of  noise)  from  a perfectly  reflecting  point  target  will  be  of  the 
form 


r(t)  = a(t-r)  cos[(w  +w.)(t-x)  + 6(t-x)]  (2.1-2) 

o a 

where  r and  correspond  to  the  target's  range  and  range  rate,  respective- 
ly. Once  the  values  of  x and  w.  are  determined,  the  point  target's  range  and 

d 

radial  range  rate  are  found  from  the  relationships 

R = ct/2  , R = coi,/2  w 

d o 

where  c denotes  the  velocity  of  light. 

Although  expression  (2.1-2)  yields  a generally  good  model  for  the  de- 
terministic component  of  the  return  signal  from  a point  target,  it  must  be 
appreciated  that  the  actual  return  signal  will  be  contaminated  by  additive 
noise.  A more  realistic  return  signal  model  will  be  then  given  by 

x(t)  = r (t)  + n(t) 

= a(t-x)  cos[(u  -hD,)(t-x)  + 8(t-x)]  + n(t) 
o a 

where  n(t)  denotes  a nondeterministic  noise  signal  which  can  arise  from  such 
sources  as  atmospheric  conditions,  clutter,  chaff,  etc.  The  main  objective 
of  radar  signal  processing  is  that  of  processing  this  noise  contaminated 
signal  so  as  to  obtain  an  appropriately  accurate  estimation  of  the  delay- 
Doppler  shift  pair  (x,md).  As  might  be  anticipated,  this  estimation  will 
be  considerably  enhanced  if  the  modulating  signals  a(t)  and  0(t)  are  judi- 
ciously chosen. 

It  will  be  beneficial  to  depict  the  radar  signal  processing  model  as 
shown  in  Figure  2.1-1.  The  transmitted  signal  s(t)  is  shown  as  being  oper- 
ated upon  to  produce  the  corresponding  point  return  signal  model  r(t)  which 
in  turn  has  added  to  it  the  nondeterministic  noise  signal  n(t).  This  sum 
signal  is  then  operated  upon  by  a radar  signal  processor  to  generate  esti- 
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s(t)  = a(t)  cos[to  t] 
o 

r(t)  - a(t-T)  cos[  (wo+a)d)  (t-T)  + 0(t-t)] 

FIGURE  2.1-1  MODEL  OF  RADAR  SIGNAL  PROCESSING 

2.2  TARGET  RESOLUTION  AND  THE  AMBIGUITY  FUNCTION 

One  of  the  primary  tasks  of  a radar  signal  processor  is  that  of  detec- 
ting the  presence  of  (or  resolving)  two  or  more  targets  in  the  return  signal. 
As  is  well  known,  this  capability  is  very  much  dependent  on  the  selection  of 
the  envelope  and  phase  modulating  signals  a(t)  and  0(t),  respectively  of  the 
transmitted  signal.  This  is  readily  made  apparent  by  considering  the  special 
case  in  which  there  are  two  point  targets  in  space  that  have  the  associated 
delay-Doppler  shift  pairs  and  (tj+T . The  two  targets  are  then 

said  to  be  separated  by  the  differential  delay  r and  differential  Doppler 
shift  to.  Using  the  previous  section's  postulated  point  target  model,  it 
follows  that  the  noise  free  return  signal  from  each  of  these  point  targets 


(2.2-lb) 


r,(t)  » a(t-T.-T)  cos [ (uj  +u>,+ti>)  (t-T,-t)  + 0(t-t  -t)] 
t-  1 011  1 


The  ability  to  resolve  these  two  point  target  return  signals  will  be 
obviously  dependent  upon  the  distinctiveness  of  their  timewise  behavior. 

There  are  a variety  of  measures  for  quantizing  this  timewise  distinctiveness. 
Undoubtably,  the  most  widely  used  such  measure  is  the  Integral  of  the  dif- 
ference squared  (integral  squared  error)  as  given  by 
00 

E(t ,w)  - J lrx(t)  - r2(t)]2dt  (2.2-2) 


This  Integral  squared  error  criterion  yields  a very  effective  means  for 
guaging  the  distinctiveness  of  two  point  target  return  signals.  In  addi- 
tion, its  utilization  leads  to  a desirable  and  tractable  analysis. 

The  modulating  signals  a(t)  and  9(t)  are  said  to  readily  resolve  the 
two  point  target  return  signals  (2.2-1)  with  differential  delay-Doppler 
shift  pair  (t,u)),  if  the  integral  squared  error  (2.2-2)  takes  on  a suit- 
ably large  positive  value.  It  is  clear  that  this  criterion  will  be  zero 
at  t=0,  (d=0  since  the  return  signals  r^(t)  and  ^(t)  are  identical  in  that 
case.  Our  basic  objectives  in  selecting  a(t)  and  0(t)  will  be  that  of 
causing  this  integral  squared  error  measure  to  be  large  positive  for  delay- 
Doppler  shift  pairs  not  at  the  origin  of  the  (t,w)  plane. 

To  evaluate  the  postulated  integral  squared  error  criterion,  one  simply 
substitutes  the  point  target  return  signal  expressions  (2.2-1)  into  rela- 
tionship (2.2-2).  Utilizing  simple  trigonometric  identities  and  the  as- 
sumption that  the  modulating  signals  a(t)  and  0(t)  are  slowly  varying  time 
signals  relative  to  the  unmodulated  signal  cos[a>ot),^  it  is  readily  shown 


This  is  commonly  refered  to  as  the  "narrow  band  signal"  assumption  and  is 
typical  of  most  radar  signals.  In  essence,  the  narrow  band  signal  assump- 
tion sets  to  zero  all  integrals  of  the  form 


x(t)  sin[w  t + 0]  dt 

V 

where  the  signal  x(t)  is  a very  slowly  varying  function  of  time  relative 
to  the  sinusoidal  signal  sin[u)Qt  +0], 
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that  the  integral  squared  error  criterion  is  given  by  (e.g.,  see  ref  [2]) 


E(t,uj) 


-J 


a2(t)dt  - |x(t,w)|  cos  [ (w^+w^+u))  t + /x(T»“)l  (2.2-3) 


The  complex  valued  function  x(T»w)  which  appears  in  this  expression  is 
known  as  the  "ambiguity  function"  which  corresponds  to  the  transmitted 
signal  and  is  formally  given  by 


x(t,w) 


h(t)p 


(t-t)e  ^wtdt 


(2.2-4) 


where  p(t)  is  the  so-called  complexed  valued  modulating  signal  specified  by 
J0(t) 


p(t)  * a(t)eJ 


(2.2-5) 


This  complex  modulating  signal  will  play  a major  role  in  much  of  what  is  to 
follow.  It  is  seen  to  be  uniquely  specified  by  the  transmitted  signal's 
modulating  envelope  and  phase  signals. 

An  inspection  of  the  above  closed-form  expression  for  reveals 

that  the  second  term  is  the  only  one  dependent  on  the  differential  delay- 
Doppler  shift  pair  (t,oj).  Furthermore,  for  virtually  all  radar  signals  of 
interest,  the  ambiguity  function  x(T»u)  will  be  a much  slower  varying  func- 
tion of  t than  will  its  multiplicative  cosine  term.  With  this  in  mind,  it 
then  follows  that  if  one  desires  E(x,w)  to  be  as  large  as  possible  for  a 
continuous  range  of  differential  delay  t values,  it  is  necessary  that  the 
quantity  |x(t»“0|  which  multiplies  the  rapidly  varying  cosine  term  must  be 
near  zero  for  all  values  of  (t,<d)  other  than  the  origin. 

Thus,  an  ideal  ambiguity  function  would  be  one  which  is  everywhere 
zero  in  the  (t,w)  plane  except  at  the  origin  where  it  has  value  equal  to 
the  energy  in  the  signal  a(t),  that  is 


X(t,<*>) 


ideal 


= < 


-J 

0 


a (t)dt 


T=(i)=0 

otherwise 
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Although  it  is  impossible  to  generate  a physically  realizable  set  of  modu- 
lating signals  a(t)  and  0(t)  which  will  attain  this  ideal  objective,  this 
behavior  will  serve  as  a comparison  reference  for  judging  the  worth  of  prac 
tical  modulating  signals  in  resolving  multiple  targets.  In  summary  a good 
modulation  signal  selection  p(t)  “ a(t)e^^C^  is  one  whose  corresponding 
ambiguity  function  closely  approximates  this  spike-like  behavior. 

The  ambiguity  function  as  defined  by  relationship  (2.2-4)  possesses  a 
number  of  interesting  analytical  properties  which  will  characterize  its  be- 
havior in  the  (t,ui)  plane.  Some  of  the  more  basic  properties  are: 

00 

Property  1 X(0,0)  - J a2(t)dt  > |x(t,u))| 

|x(t,w) |2dTdw  = x2(0,0) 

Property  3 |x(f,u>)|  = |x(-t,-ui)| 

Property  2 is  particularly  noteworthy  in  that  it  indicates  that  the  postu- 
lated ideal  ambiguity  function  is  not  achievable.  Namely,  since  the  vol- 

2 

ume  under  the  |x(t,w)|  surface  is  a constant,  any  attempt  to  cause  |x(t,u>) 
to  be  spike-like  will  invariably  result  in  a large  spread  of  nonzero  levels 
off  the  origin. 

2.3  THE  MATCHED  FILTER 

One  of  the  primary  objectives  of  any  radar  signal  processor  is  that  of 
detecting  the  existence  of  a weak  return  signal  in  the  presence  of  additive 
noise.  This  detection  will  be  enhanced  if  the  linear  processor's  unit- 
impulse  response  is  chosen  to  in  some  sense  match  the  return  signals  wave- 
form. To  illustrate  this  concept,  let  us  consider  the  specific  case  in 
which  the  input  is  given  by 

x(t)  - r(t)  + n(t) 

where  r(t)  denotes  the  useful  return  signal  component  and  n(t)  is  additive 


Property  2 


1_ 

2tt 


r 
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noise.  If  the  additive  noise  is  white,  it  is  well-known  that  the  "matched 
filter"  whose  unit-impulse  response  is  specified  by 

h(t)  = r*(-t)  (2.3-1) 

will  maximize  the  output  signal -to-ioise  power  ratio  at  time  t=0.^  A linear 
filter  with  this  unit-impulse  response  is  then  said  to  be  optimally  matched 
for  detecting  the  signal  r(t)  which  is  masked  by  additive  white  noise.  Sim- 
ilar results  may  be  obtained  when  the  additive  noise  is  other  than  white 
(e.g. , see  ref.  [2])  . 

With  this  in  mind,  let  us  consider  the  typical  matched  filter  receiver 
shown  in  Figure  2.3-1  in  which  the  input  signal  is  taken  to  be  of  the  form 

r(t)  = a(t-t)  cos[  (o)o+w)  (t-r)  + 9(t-t)]  (2.3-2) 

and  where  the  delay-Doppler  shift  pair  (t,w)  is  not  known  apriori.  It  will 
be  assumed  that  the  matched  filter  was  designed  to  detect  the  presence  of  a 
return  signal  whose  Doppler  shift  is  io^.  The  unit-impulse  response  of  the 
required  matched  filter  is  readily  found  to  be 

h(t)  = a(-t)  cos[(w  -w  -ho, ) t-0 (-t)J  (2.3-3) 

o a i 

It  can  be  shown  that  the  response  of  this  matched  filter  to  input  signal 
(2.3-2)  is  given  by 

yc(t)  = | I x(T-t,u)-u)1)  I cos[(u)o-u)a+u))t  + / x(Tt,oi-m1)  ] 
where  x(t5m)  is  the  ambiguity  function  corresponding  to  the  complexed  valued 

^Generally,  for  purposes  of  physical  realizability,  the  required  matched 

filter's  unit-impulse  response  must  be  appropriately  time  shifted  for  causal 

* 

signal  processing  (i.e.,  r (T^-t)).  In  this  case,  the  output  signal  to  noise 
power  ratio  is  maximized  at  t=T  . 
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modulating  signal  p(t)  = a(t)e^e^.  In  arriving  at  this  result,  the 
narrow-band  signal  assumption  is  made  in  which  the  intermediate  frequency 
is  taken  to  be  very  large. 

Thus,  the  signal  processor's  output  is  seen  to  have  an  envelope  which 
is  equal  to  a translated  version  of  the  transmitted  signal's  associated 
ambiguity  function  magnitude.  If  the  matched  filter  is  "tuned"  to  the  in- 
put signal  (i.e.,  u>  * u) , then  the  envelope  of  the  output  signal  will  peak 
at  time  t=T  thereby  indicating  the  presence  of  a target  at  the  range  corres- 
ponding to  delay  t.  This  particular  form  of  matched  signal  processing  has 
served  as  a very  effective  means  for  detecting  the  presence  of  weak  target 
return  signals  in  a background  of  strong  noise.  Namely,  one  looks  for  the 
presence  of  pulse  like  peaks  in  the  envelope  and  measures  the  times  at  which 
they  occur.  From  the  above  remarks,  it  is  apparent  that  the  ambiguity  func- 
tion takes  on  added  meaning  whenever  matched  filtering  operations  are  uti- 
lized. 


FIGURE  2.3-1  MATCHED  FILTER  RECEIVER 


2.4  AUTOCORRELATION  FUNCTION 

The  ability  to  achieve  fine  range  measurements  is  inherently  linked 
to  the  time  behavior  of  the  autocorrelation  function  which  is  associated 
with  the  modulating  signal.  This  is  readily  demonstrated  by  noting  that 
the  envelope  of  the  response  of  matched  filter  (2.3-3)  to  the  input  signal 
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given  by  expression  (2.3-2)  with  w - is  simply 


ye(t)  = 4 lx(T-t,  0) 


(2 . 4-1) 


The  ambiguity  function  evaluated  at  u=0  is  otherwise  known  as  the  auto- 
correlation function  associated  with  p(t)  and  is  formally  defined  by 


^(t)  = x(t,0) 


■00 

li(t)u 

o< 


(t-T)dt 


(2.4-2) 


The  value  of  this  autocorrelation  function  at  t=0  is  seen  to  be  equal  to  the 

energy  contained  in  the  modulating  signal.  It  is  a well-known  fact  that  the 

autocorrelation  function's  magnitude  is  maximum  at  t=0  (i.e.,  R (0)  > |r  (t) I 

P — p 

for  all  r). 

In  summary,  we  have  established  the  following  very  important  input- 
response  pair  relationship  for  the  matched  filter  operator  (2.3-3) 


x ( t ) = a(t-t)  cos[ (u^+u^) (t— t)  + 9 ( t— t ) ] 


yjt)  = t R (t-T) 
e 4 p 


(2.4-3) 


To  appreciate  the  significance  of  this  relationship,  let  us  consider  those 
modulating  signals  p(t)  which  have  autocorrelation  functions  which  are  every- 
where almost  zero  except  near  the  origin  (i.e.,  spike-like  in  appearance).  For 
such  modulating  signals,  the  matched  filter  response  (2.4-3)  is  a signal 
which  is  everywhere  zero  except  near  t=T  where  a spike  will  appear.  Thus, 
to  detect  the  presence  of  target  echoes  in  the  return  signal,  one  would 
simply  examine  the  matched  filter's  response  signal  for  the  existence  of 
spikes.  The  existence  of  a spike  would  indicate  the  presence  of  a target 
whose  range  is  proportional  to  the  time  occurance  of  that  spike.  Clearly, 
the  narrower  the  width  of  the  spike,  the  more  precisely  the  target's  range 
can  be  determined. 

With  the  above  thoughts  in  mind,  it  is  apparent  that  the  class  of  modu- 
lating signals  which  have  impulse-like  autocorrelation  functions  are  of  great 
usefulness  in  radar  signal  processing.  The  ideal  autocorrelation  function 
would  be  one  for  which  R^(t)  is  identically  zero  except  at  t=0  where  its 


value  equals  the  energy  In  the  modulating  signal.  Although  this  ideal  ob- 
jective cannot  be  achieved  by  physically  realizable  modulating  signals,  it 
serves  as  a basis  for  comparing  the  effectiveness  of  different  modulating 
signals. 


SECTION  3 


PULSE -BURST  SIGNALS 


3.1  BASIC  DEFINITIONS 

A particularly  effective  method  for  generating  radar  signals  which 
possess  desirable  fine  range  resolution  characteristics  is  that  of  pulsed 
sinusoidal  modulation.  In  general,  a pulsed  sinusoid  consists  of  a clust- 
ered sequence  of  L+l  sinusoid  signals  each  of  duration  T seconds.  In  each 
of  these  L+l  "inner-pulse"  intervals,  the  sinusoid's  amplitude  and  phase 
angle  will  be  held  fixed,  but,  these  parameters  will  generally  vary  from 
one  inner-pulse  interval  to  the  next.  Signals  so  characterized  are  commonly 
referred  to  as  "pulse -burst  signals"  and  have  the  following  convenient  ana- 
lytical representation 

L 

s ( t ) = 1 a cos [m  t+0  ] (u(t-kT)-u(t-kT-T) } (3.1-1) 

k=0  K 0 K 

where  a^  and  6^  designate  the  transmitted  sinusoid's  amplitude  and  phase 
angle,  respectively,  during  the  kth  inner  pulse  interval  kT  _<  t < kT+T  and 
u(t)  denotes  the  standard  unit-step  signal.''  This  pulse-burst  signal  is 
seen  to  have  a time  duration  of  LT+T  seconds  and  is  identically  zero  out- 
side the  time  interval  0 _<  t < LT+T.  We  shall  refer  to  this  pulse  burst 
signal  as  having  length  L+l  and  its  time  behavior  is  seen  to  be  completely 
specified  by  the  2L+2  real  numbers  (a^,  a^,  . . . , a^}  and  {0^,  0^,  . . . ,0^}. 

In  order  to  determine  the  suitability  of  this  signal  for  radar  proc- 
essing applications,  let  us  first  obtain  the  complex  modulating  signal  p(t) 
which  corresponds  to  it.  According  to  expressions  (2.1-1)  and  (2.2-5),  this 
complex  modulating  signal  is  specified  by  the  following  piecewise  constant 
function  of  time 


1 


The  unit-step  signal  is  defined  by 


u(t) 


t > 0 
t < 0 
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(3.1-2) 


p(t)  = l a e k(u(t-kT)  - u(t-kT-T)} 
k-0  k 

Again,  it  is  clear  that  magnitude  and  phase  angle  parameters  {a^}  and  { 0 ^ } 
play  a major  role  in  characterizing  the  underlying  radar  signal.  In  par- 
ticular, the  associated  ambiguity  function  (2.2-4)  which  depicts  the  de- 
sirabiltiy  of  the  transmitted  signal  for  radar  processing  is  seen  to  be 
explicitly  dependent  on  y(t)  and  therefore  on  these  parameters. 

PULSE- BURST  CODE 

The  suitability  of  a given  pulse-burst  signal  for  radar  processing  is 
then  completely  determined  by  the  selection  of  the  L+l  complex  numbers 

j°k 

uk  = ake  for  k = 0,  1,  ....  L (3.1-3) 

It  is  clear  that  a knowledge  of  these  elements  uniquely  identifies  the 
pulse-burst  signal  (3.1-1)  and  vice-versa.  As  such,  the  finite  sequence 

{uo’  ui’  • • * » (3.1-4) 

will  be  refered  to  as  the  "pulse  burst  code"  associated  with  the  pulse 
burst  signal  (3.1-1).  In  much  of  what  is  to  follow,  we  will  be  concerned 
with  presenting  (or  developing)  methods  for  selecting  the  code  elements 
so  that  the  autocorrelation  function  associated  with  the  complex  modu- 
lating signal  p(t)  takes  on  desirable  characteristics. 

AMBIGUITY  FUNCTION  OF  A PULSE-BURST  SIGNAL 
It  has  been  previously  indicated  that  the  ambiguity  function  associated 
with  a given  complex  modulating  signal  provides  an  attractive  visual  method 
for  judging  a signal's  effectiveness  in  radar  processing.  With  this  in  mind, 
the  ambiguity  function  which  corresponds  to  pulse  burst  signal  (3.1-1)  was 
next  found  and  is  given  by 
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• • • M 


X(nT+A,w)  = (T-A)e 


(3.1-5) 


1u»(T-A) 


L 


2 , roj(T-A)  , v * -jwkT  . 

sine  l—  1 l Vt_„e 

k=n 


ju)(2T-A) 


Ae 


Sinc[2"]  . I . Vk-n-1 

k-n-f  1 


* e-jcokT 


for  n = 0,  ±1,  ±2,  ±3,  . 


and  0 < A < T 


An  examination  of  this  ambiguity  expression  indicates  two  interesting  char- 
acteristics of  pulse-burst  signals.  First  of  all,  this  ambiguity  function 
is  identically  zero  for  delays  greater  than  the  pulse-burst  length  LT+T 
(i.e.,  | n | > L) . Furthermore,  since  the  two  summations  in  this  relationship 
are  each  periodic  functions  of  w with  period  2ir/T,  it  follows  that  the  am- 
biguity function  will  tend  to  periodically  decay  to  zero  in  a sin  x/x  manner 
in  the  w direction  for  a fixed  delay  time. 

Once  the  values  for  the  pulse  code  p^,  p , . . . , have  been  assigned, 
it  is  a simple  matter  to  generate  the  surface  behavior  of  the  ambiguity  func- 
tion (3.1-5)  with  the  aid  of  a suitable  computer  program  and  plotting  routine 
(e.g.,  see  Appendix  8.1).  In  addition,  with  the  use  of  this  convenient  ex- 
pression, the  distinct  possibility  exists  that  one  may  generate  an  algorithmic 
procedure  for  selecting  the  code  elements  p^  so  that  the  ambiguity  function 
will  be  optimum  in  some  sense  (e.g.,  impulse-like  in  behavior).  This  will 
be  demonstrated,  to  some  extent,  in  the  next  section  when  these  code  elements 
are  chosen  to  make  the  associated  autocorrelation  function  impulse-like. 


3.2  AUTOCORRELATION  OF  A PULSE-BURST  SIGNAL 

As  indicated  in  Section  2. A,  an  important  objective  in  many  radar  signal 
processing  applications  is  the  generation  of  finite  duration  modulating 
signals  whose  autocorrelation  is  impulse -like.  The  autocorrelation  function 
which  corresponds  to  a pulse  burst  signal  is  readily  obtained  by  evaluating 
the  associated  ambiguity  function  (3.1-5)  at  m=0  to  yield 

i 

1 

l 

1A 


l 


(3.2-n 


R (nT+A) 
M 


(T-A) 


L 

l 


„ Vk-n  + * . I . VWl 

k=*n  k-n+l 


n=0,  ±1,  ±2,  . . . and  0 < A < T 


From  this  expression,  it  is  apparent  that  this  autocorrelation  function 
R (t)  is  identically  zero  outside  the  interval  |t|  < LT+T.  This  is  an 
obvious  consequence  of  the  fact  that  the  duration  of  the  modulating  signal 
p(t)  is  LT+T  seconds. 

It  will  be  desirable  to  investigate  this  autocorrelation  function  at 
integral  multiple  values  of  T seconds  (i.e.,  A=0) . This  yields  the  sampled 
autocorrelation  sequence 


R^(nT)  = T 


L 

l 

k=n 


Vk-n 


(3.2-2) 


whose  elements  are  identically  zero  for  integer  values  n satisfying  |n|  > L. 
It  is  possible  to  reconstruct  the  original  autocorrelation  function  (3.2-1) 
from  this  associated  sampled  autocorrelation  function  by  means  of  the  simple 
relationship 


R (nT+A)  = (l-A)R  (nT)  + AR(nT+T)  for  0 < A < T 

M V ~ 

This  is  observed  to  be  a straight  line  interpolation  between  the  sampled 
values.  Thus,  pulse-burst  signals  as  represented  by  modulating  signal 
(3.1-2)  will  have  associated  autocorrelation  functions  which  are  piecewise 
linear  and  continuous. 


3.3  QUALITY  MEASURES  OF  A PULSE -BURST  SIGNAL 


The  code  sequence  p^,  p , 


p which  corresponds  to  a pulse-burst 
Li 


signal  is  said  to  possess  an  impulse-like  autocorrelation  if  the  sampled 
auotcorrelation  is  such  that  the  normalized  maximum  side  lobe  as  defined  by 


A = max 
n^O 


|Ry(nT) 

R (0) 
h 


(3.3-1) 
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is  sufficiently  close  to  zero  (note,  0 <•  A _<  1).  In  other  words,  the  zero''1' 
element  of  an  impulse-like  sampled  autocorrelation  sequence  dominates  the 
magnitude  of  its  remaining  elements.  A variety  of  systematic  procedures 

for  selecting  the  code  sequence  pQ,  u have  been  developed  for 

achieving  such  autocorrelation  function  behavior.  Before  presenting  two 
of  the  more  popular  of  such  methods,  let  us  first  outline  some  practical 
considerations  which  are  imposed  on  pulse-burst  signals. 

In  order  to  achieve  satisfactory  probability  of  detection  for  a given 
false  alarm  rate,  it  is  necessary  that  sufficient  signal  energy  be  trans- 
mitted. For  pulse-burst  signals,  the  transmitted  signal's  energy  is  equal 
to  one  half  the  energy  in  the  modulating  signal  p(t).  The  modulated  signals 
energy,  however,  is  simply  equal  to  R(0),  that  is 


Vk 


(3.3-2) 


This  quantity  must  be  made  larger  than  some  prescribed  value  in  order  to 
obtain  a satisfactory  probility  of  detection  for  a given  false  alarm  rate. 
There  are  three  methods  for  increasing  the  modulating  signal's  energy 


(i)  increase  the  inner- pulse  duration  time  T 

(ii)  increase  the  magnitude  of  the  code  elements  p^ 

(iii)  increase  the  length  of  the  pulse  code 


Since  the  main  spike  of  the  pulse-burst  autocorrelation  function  is  defined 
over  the  time  interval  -T  _<  t _<  T,  fine  range  resolution  considerations  re- 
strict the  size  of  the  inner-pulse  duration  time  T.  Furthermore,  peak  power 
limitations  of  the  radar  transmitter  restrict  the  magnitude  of  the  individual 
code  elements  to  be  no  larger  than  a fixed  quantity.  With  this  in  mind,  it 
follows  that  the  most  practical  method  for  increasing  the  energy  of  the 
transmitted  signal  is  to  increase  the  code  length  and  to  simultaneously 
have  each  of  the  code  element  values  p^  be  nearly  equal  in  magnitude  (i.e., 
peak  power  is  nearly  maintained  throughout  the  burst). 


In  selecting  the  elements  p^,  p , 


p^  of  a pulse-burst  code,  wc 
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will  then  be  concerned  with  both  obtaining  an  autocorrelat ion  function  which 
is  impulse-like  while  simultaneously  maintaining  a reasonably  efficient  means 
of  transmitting  signal  energy.  A criterion  which  measures  this  last  criterion 
is  the  so-called  energy  efficiency  as  defined  by  [5] 


P 


burst  signal  energy  transmitted 

energy  transmitted  at  peak  power  over  each  T second  duration 


(3.3-3) 


(L+l)  (max | pj  2) 


E 

y 

(L+l) (max| PjJ  2) 


This  criterion  takes  on  values  in  the  interval  0 ■'  p j 1 with  the  most 
energy  efficient  codes  being  those  with  equal  code  magnitude  and  p=l.  To 
measure  a codes  variance  from  the  most  efficient  energy  mode  (i.e.,  p=l), 
the  following  criterion  lias  been  also  used  [51 


e 


(3.3-4) 


Although  a given  pulse-burst  code  may  possess  an  autocorrelation  function 
which  is  impulse- like,  it  may  be  deemed  unsatisfactory  if  its  energy  effi- 
ciency is  too  small.  A truly  good  pulse-burst  code  is  one  whose  individual 
elements  are  nearly  equal  in  magnitude  and  which  has  an  impulse -like  auto- 
correlation function.  The  Barker  codes  to  be  presented  in  the  next  section 
are  an  example  of  such  a code. 


3.4  CONSTANT  MAGNITUDE  CODES 

When  transmitting  a pulsed -sinusoidal  signal  under  amplitude  constraints, 
maximum  power  transmittal  is  achieved  by  maintaining  a maximum  sinusoid 
amplitude.  In  terms  of  a pulse-burst  code,  this  corresponds  to  having  t lit* 
elements  of  the  code  be  specified  by 

j9k 

Pk  = Ae  k = 0,  1,  . . . , L (3.4-1) 
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where  A corresponds  to  the  peak  amplitude  limitation.  The  energy  utilization 
efficiency  for  a code  of  this  nature  is  optimum  (i.e.,  p=l)  and  its  variance 
is  likewise  minimal  (i.e.,  c=0) . We  shall  refer  to  a code  whose  elements  are 
speeit  ied  hy  relationship  (3.4-1)  as  a constant  magnitude  code  to  ret  loot  t lie 
fact  that  |a^j=A  for  k = 0,  1,  2,  . . . ,1.. 

In  synthesizing  a constant  magnitude  code,  one  seeks  to  select  the  phase 
terras  0^,  0^,  . . . , 0^  so  that  the  resultant  transmitted  signal  has  an  as- 
sociated autocorrelation  function  which  is  near  impulse-like  in  behavior. 

The  autocorrelation  sequence  which  corresponds  to  code  (3.4-1)  is  readily 
obtained  by  utilzation  of  expression  (3.2-2)  and  results  in 


R (nT)  = AT 
P 


L 

I 

k=n 


JiVVJ 


n = 0,  1, 


(3.4-2) 


BINARY  CODES  AND  THE  BARKER  CODES 


One  of  the  most  widely  studied  of  constant  magnitude  codes  are  those 
which  the  phase  angle  elements  8^  take  on  the  exclusive  values  of  0 or 
Such  codes  are  also  called  binary  codes  in  that  the.  code  elements  p 


take  on  either  of  the  two  real  values  A or  -A.  Binary  codes  are  of  inter- 
est mainly  because  of  their  optimum  energy  utilization  efficiency,  and,  the 
ease  with  which  they  may  be  generated. 

To  illustrate  a particular  important  class  of  binary  codes,  lot  us  now 
consider  the  Barker  codes.  The  code  (p^,  P ^ , p^)  is  said  to  be 
Barker  if  its  elements  are  exclusively  plus  or  minus  one  and  its  associated 
autocorrelation  sequence  is  specified  by 

r L+l  n = 0 

(3.4-3) 

0 or  +1 


R (nT)  = i 

P 


n = 0 
n ^ 0 


It  is  observed  that  the  normalized  maximum  side  lobe  ratio  for  a Barker  code 
is  1/(L+1) . Thus,  Barker  codes  would  be  ideal  vehicles  for  obtaining  impulse- 
like  autocorrelation  functions  if  the  pulse-burst  length  L+l  could  be  made' 
suitably  large.  Unfortunately,  Barker  codes  exist  only  for  a small  set  of 
lengths  as  shown  in  Table  3.4-1.  Thus,  Barker  codes  have  somewhat  limited 
use  and  are  generally  not  particularly  appropriate  in  applications  requiring 


near  impulse-like  autocorrelation  belavior. 


Code  Length 


TABLE  3.4-1 


Code  Elements 
1 

1,  -1 
1.  1.  “I 
1.  1,  1,  -1 

1.  1.  1,  "I,  1 

1,  1,  1,  -1,  -1,  1,  -1 

1,  1,  1,  -1,  -1,-1,  1,  -1,  -1,  1,  -1 

1,  1,  1,  1,  1,  -1,  -1,  1,  1,  -1,  1,  -1,  1 

BARKER  CODES 


To  illustrate  the  characteristics  of  a Barker  code,  let  us  consider  the 
length  11  code  as  depicted  in  Figure  3.4-1.  The  autocorrelation  function's 
magnitude  which  corresponds  to  this  code  is  also  shown  along  with  the 
parameters  characterizing  the  code's  efficiency.  In  a later  section,  we 
will  compare  this  code's  characteristics  with  those  of  a so-called  pseudo- 
optimum  code  generated  by  an  algorithm  to  be  presented  in  Section  4. 
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3.5  HUFFMAN  CODES 


In  generating  pulse-burst  codes,  superior  lmpulse-like  autocorrelation 
behavior  can  be  achieved  by  allowing  the  code  elements  p^  to  be  uncontrained 
complex  numbers.  It  is  of  course  realized  that  in  utilizing  such  codes, 
particular  attention  has  to  be  paid  to  their  associated  energy  utilization 
efficiency  factor  p when  peak  power  limitations  are  imposed.  The  class  of 
Huffman  codes  to  be  now  considered  offer  an  excellent  example  of  high  per- 
formance pulse-burst  codes. 

The  generally  complex  valued  code  {pQ,  p^,  . . . , u^}  is  said  to  be 
a Huffman  code  if  its  associated  autocorrelation  sequence  is  given  by 


V”T> 


n » 0 
n = ±L 
otherwise 


(3.5-1) 


where  = E|pjJ  denotes  the  energy  in  the  code  sequence.  This  auto- 
correlation sequence  can  be  made  more  impulse-like  by  simply  increasing 
the  code's  energy  E^.  As  a further  observation,  it  is  noted  that  the  re- 
quirement R (LT)  “ p p - -1  imposes  a constraint  on  the  code's  first  and 

|i  U Li 

last  element  values. 

In  what  is  to  follow,  we  shall  be  concerned  with  presenting  a systematic 
procedure  for  selecting  a L+l  length  complex  code  which  will  have  the  specfied 
autocorrelation  behavior  (3.5-1).  This  is  best  begun  by  considering  the 
z-transform  of  the  finite  length  code  {p^}  as  defined  by 


-1  -2 

P(z)  = Pq+P-^z  +P2Z  + • • • +\z 


(3.5-2) 


where  z is  a complex  variable.  Noting  that  P(z)  is  a polynomial  in  the 
variable  z \ it  follows  that  the  code  transform  has  the  equivalent  factorej 
representation 


*The  Huffman  code  in  which  R (LT)  =»  e^6 

P 

similar  to  that  to  be  presented  in  this 


can  be  synthesized  in  a manner 
section. 
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A 


(3.5-3) 


P(z)  - pQz  L(z-z1) (z-z2)  . . . (z-zL) 


where  z,  are  the  roots  of  P(z) 
k 

in  some  manner,  the  set  of  roots  z^,  z^, 


Clearly,  if  it  is  possible  to  determine. 


z^  of  P(z)  and  the  con- 


stant Uq,  then  one  may  generate  the  corresponding  code  elements  by  carrying 
out  the  indicated  multiplications  in  expression  (3.5-3).  Such  a procedure 
shall  now  be  presented. 

The  z-transform  of  the  autocorrelation  sequence  corresponding  to  the 


code  {uQ,  p^ 


S(z) 


, p^)  is  readily  found  to  be 


L 

L „ 

c -k 

p * n 

l 

_ k«0  K 

V 

n*0 

(3.5-4) 


-L 


L 

* 

[ L * I 

n 

k-l 

J 

i 

n (z-i/z  ) 

i n 

n°l 

-z 


where  use  of  the  constraint  p p.  ■ -1  has  been  made.  It  is  seen  that  the 

U L 

roots  of  S(z)  occur  in  reciprocal  complex  conjugate  pairs  with  L of  the 
roots  belonging  to  P(z)  and  their  corresponding  L reciprocal  complex  con- 
jugate pair  mates  belong  to  the  polynomial  multiplying  P(z)  in  expression 
(3.5-4). 

If  a particular  complex  code  sequence  {p^,  p^,  ....  p^}  is  to  hav 
the  specified  autocorrelation  sequence  (3.5-1),  these  code  elements  must 
then  be  selected  such  that 


S(z) 


-L.  2L  _ L . 
-z  (z  -E  z +1) 


(3.5-5) 


2L  L 

The  polynomial  z -E^z  +1»  which  in  part  constitutes  the  transform  S(z),  is 
observed  to  be  a mirror- image  polynomial  with  real  coefficients.  It  follows 


Id  16 

that  if  reJ  is  a root  of  S(z),  then  so  must  reJ  , 


r ^e^0,  and,  r ^ ^ also 


where  the  scalar  r is  given  by 


Thus,  these  2^  roots  are  seen  to  lie  on  two  concentric  circles  centered  at 
the  origin  along  rays  spaced  every  2ir/L  radians.  This  root  pattern  is 
depicted  in  Figure  3.5-1  for  the  special  case  N-8. 


To  generate  the  roots  of  polynomial  P(z),  one  then  simply  selects  z^ 
to  be  equal  to  either  re^fflc^L  or  r for  ^ - 1,  2,  ....  L. 

Carrying  out  this  sequential  root  selection  results  in  the  following  code 
transform 

P(z)  = az  L(z-z1) (z-z2)  . . . (z-zL)  (3.5-7) 

It  is  observed  that  the  L roots  not  selected  In  this  manner  correspond  to 

the  reciprocal  complex  conjugate-pair  mates  of  the  set  {z.,z  , . . . , z }. 

I z * L 

The  constant  a is  chosen  so  as  to  satisfy  the  imposed  condition  p p = -1, 

U L 

namely 

1 L 

1*1  n (“*k>  " “I 

k-1  K 

As  a final  step,  the  multiplications  in  expression  (3.5-7)  are  carried  out 
to  yield  the  code  element  values  {p^}- 

From  this  construction  process,  it  is  apparent  that  since  there  are 
two  choices  for  each  of  the  L roots  of  P(z)  on  each  half  ray,  there  will 
then  exists  a total  of  2L  different  code  sequences  which  will  have  the 
same  autocorrelation  sequence  (3.5-1).  As  indicated  in  reference  [5], 
some  of  these  seemingly  different  codes  are  in  fact  equivalent  in  terms  of 
the  code  magnitude  distribution  (i.e.,  |p^|,  |p2I,  • * • » I | ) - More- 
over, some  of  the  code  selections  will  have  a relatively  good  energy  ef- 
ficiency p while  others  will  be  very  unsatisfactory  in  this  regard.  In 
addition,  the  ambiguity  functions  which  correspond  to  the  individual  2L 
codes  will  be  highly  variable  in  terms  of  desired  spike-like  behavior  in 
the  (t,<d)  plane.  Unfortunately,  for  large  values  of  code  length  L+l,  an 
investigation  of  each  of  the  2L  realizations  is  not  feasible.  Some  rel- 
atively ad-hoc  methods  based  on  conjectures  have  been  offered  for  obtaining 
individual  satisfactory  realizations.  For  example,  the  "optimum"  length 
11  code  in  terms  of  energy  efficiency  is  depicted  in  Figure  3.5-2  along 
with  its  associated  autocorrelation  function  [5]. 

Before  leaving  the  interesting  subject  of  Huffman  codes,  let  us 
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consider  the  subset  of  the  2 selections  which  will  result  in  real  codes. 

A little  thought  will  indicate  that  real  codes  will  result  if  and  only  if 

the  roots  z constituting  P(z)  are  selected  to  occur  in  complex  conjugate" 

^ * 

pairs  (i.e.,  z = z ).  When  L is  even,  there  will  exist  2 real  roots 

L~fc  * (L-2) /2 

and  L-2  complex  roots  thereby  indicating  a total  of  2 root  selections 

will  result  in  real  codes.  Similarly,  a total  of  2^L  different  real 

codes  can  be  generated  when  L is  odd. 


SECTION  4 


I 


ALGORITHMIC  GENERATION  OF  REAL  PULSE-BURST  CODES 


4.1  INTRODUCTION 

In  this  section,  we  shall  direct  our  attention  exclusively  to  bi- 
phase codes  in  which  the  phase  elements  are  either  0 or  it . A bi-phase 
code  is  typically  referred  to  as  a real  burst  code  due  to  the  fact  that 
the  associated  complex  modulating  sequence  . . . ,PL)  consists 

entirely  of  real  numbers.  It  is  to  be  noted  that  the  Barker  codes  are 
real  while  a subset  of  Huffman  codes  are  themselves  real.  Interest  in  real 
pulse*burst  codes  arises  due  to  the  simplicity  of  their  generation  as  a 
modulated  sinusoid  since  the  only  variable  is  the  sinusoid's  amplitude  and 
not  its  phase. 


4.2  PSEUDO-OPTIMUM  CODE 

A pulse-burst  code  Is  said  to  be  ideal  if  its  associated  sampled 
auotcorr elation  sequence  is  equal  to 


Ri(n)  - E6(n) 


n=0 

n?*0 


(4.2-1) 


in  which  6(n)  designates  the  unit-Kronecker  sequence  and  E denotes  the 
code's  energy^.  It  is  not  difficult  to  show  that  the  only  finite  length 
codes  which  will  have  this  characteristic  are  of  length  one.  Unfortunately, 
length  one  codes  are  not  of  Interest  since  the  main  purpose  for  using 
pulse*burst  codes  is  that  of  transmitting  a long  duration  pulse  (to  over- 
come peak  power  limitations)  which  in  turn  can  be  effectively  transformed 
into  a short  duration  pulse  by  means  of  a matched  filtering  operation. 
Although  the  above  ideal  autocorrelation  behavior  cannot  be  achieved  by 


1 


The  unit 


«(n) 


Kronecker  sequence  is  defined  by 

0 n-0 

= « 

1 n-±l,±2,i3,  . . . 
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a finite  code  of  length  greater  than  one,  it  can  be  closely  approximated 
if  the  code  elements  are  judiciously  chosen.  We  shall  now  direct  our 
attention  to  this  most  relevant  task. 

In  what  is  to  follow,  we  shall  investigate  the  sampled  autocorrelation 
sequence  which  corresponds  to  the  L+l  length  real  pulse-burst  code 
{ V»o » • • • » In  referring  to  the  code  as  having  length  L+l,  it 

is  tacitly  implied  that  its  two  end  elements  p^  and  p^  need  to  be  each 
nonzero  for  otherwise  the  code  would  be  of  a shorter  length.  To  measure 
how  closely  this  code's  sampled  autocorrelation  sequence  R^(n)  approximates 
the  ideal  objective  (4.2-1),  the  standard  squared  error  criterion  will  be 
used,  that  is 

L , 

f(p)  - l w(n) (R  (n)  - R (n) T (4.2-2) 

n— L y 

9 L 9 

- w(0){R  (O)-EJ  + 2 l w(n)R  (n) 


The  elements  w(n)  are  nonnegative  numbers  which  are  used  to  weight  the 
approximation  errors  in  any  desired  manner.  Due  to  the  symmetry  of  the 
autocorrelation  sequence,  no  loss  of  generality  is  introduced  by  assuming 
that  w(-n)  * w(n) . As  a final  observation,  it  is  apparent  from  the  formula 
for  R^(n)  as  given  by 


Rp(n) 


l Vk-n 

k*n 


that  the  above  criterion  f(p)  is  a quartic  function  of  the  code  elements. 

It  is  now  desired  to  construct  an  acceptable  code  which  will  make  this 
squared  error  criterion  minimum.  Unless  a constraint  of  some  form  is 
imposed  on  this  minimization,  however,  it  is  seen  that  an  optimum  code  of 
length  one  will  cause  f(p)  to  take  on  the  value  zero  (e.g.,  Pq  ■=  fH, 
p » p0  =...=■  p ■ 0) . In  order  to  avoid  this  possibility,  we  shall 
constrain  the  first  code  elements  p^  to  take  on  a fixed  value  a and  then 


28 


select  the  remaining  L code  elements  to  minimize  criterion  (4.2-2).^  This 
constrained  minimization  problem  is  then  expressed  as 


min 


U0 


“ 

V 2 

w(0)  - 

l - E 

. k-0 

(4.2-3) 


where  the  analytic  expression  for  R^(n)  has  been  incorporated.  A code 
which  renders  this  criterion  a relative  minimum  value  will  be  referred 
to  as  a "pseudo-optimum"  code. 

In  order  to  encourage  a uniform  distribution  of  code  power  (i.e., 
each  element  has  magnitude  E/(L+1)),  a logical  selection  for  the  constant 
a is  E/(L+1).  It  will  be  shown  in  section  5,  however,  that  other  se- 
lections will  lead  to  superior  codes.  With  this  in  mind,  we  shall  let  Uq 
take  on  various  values  and  find  the  minimizing  codes  which  correspond  to 
each  of  these  \1q  selections.  From  this  family  of  psuedo-optimum  codes, 
we  shall  select  one  which  is  deemed  to  be  superior  (e.g.,  the  one  with 
maximum  energy  utilization  efficiency). 

A careful  examination  of  criterion  (4.2-3)  indicates  it  to  be  a re- 
latively highly  nonlinear  function  of  the  code  elements  u^.  A closed-form 
solution  for  a minimizing  code  is  then  not  feasible  and  one  must  appeal  to 
an  algorithmic  solution  approach.  In  the  next  two  sections,  a particularly 
effective  algorithmic  solution  procedure  will  be  presented. 


Another  method  for  achieving  this  objective  is  that  of  adding  the  functional 


8(n>  - I W - f^T> 


k“0 


to  f(p)  and  then  minimizing  f(p)  + Xg(p) . 


HUOCL 


- >v.  It  will  be  beneficial  to  interpret  the  finite 

•4'  • ...  ..-qiieiuf  R^(n)  a8  being  the  unit-impulse  response 

• - rftc-: tme  system  conf lgurat  ion . To  verify  this 
.!»•:  .u  p let  us  consider  the  subsystem  characterized  by 

» . .ar  4um1  nonrecurslve  operator 

■ a ♦ ..  a(a-l)  ♦ . . . + ^x(n-L)  (4.3-1) 

-ners  * a • n.  denote  the  input  and  response  sequences,  respectively, 

and  the  arc  onstants  which  will  be  subsequently  Identified  with  the 

element*  t a real  pulae-burst  code.  The  transfer  function  which  corres- 
ponds to  this  operator  is  then  given  by 

H(z)  ■ • • • +uL*  L (4.3-2) 

For  reasons  which  will  be  made  clear  shortly,  let  us  connect  in  cascade 
with  operator  (4.3-1)  its  t ime-anticausal  dual.  This  anticausal  dual  will 
be  characterized  by  the  nonrecursive  operator  relationship 

y(n)  - Uq  s(n)  + y1s(n+l)  + . . . + yLs(n+L)  (4.3-3) 

in  which  the  input  sequence  s (n)  is  equal  to  the  response  of  time-causal 
operator  (4.3-1).  This  cascaded  configuration  is  depicted  in  Figure  4.3-1 
where  it  is  noted  that  the  transfer  function  of  t ime-anticausal  operator 
(4.3-3)  is  simply  H(z  ^) . Using  standard  z-transform  theory  properties, 
it  follows  that  the  z-transform  of  this  configuration's  response  to  the 
Kronecker  unit- impulse  is  given  by 


Y(z)  = H(z)H(z  l) 


A comparison  of  this  expression  with  relationship  (3.5-4),  however,  indicates 
that  Y(z)  is  simply  the  z-transform  of  the  autocorrelation  sequence  R^(n). 

In  other  words,  the  impulse  response  of  the  nonrecursive  operator  con- 
figuration shown  in  Figure  4.3-1  is  equal  to  the  autocorrelation  sequence 
which  corresponds  to  the  real  pulse-burst  code{y  ,y  , . . . ,y  }. 

U X L 


x(n) 


H(z)  - 1 + v2z  2 + • • • + uLz  L 

FIGURE  4.3-1  CASCADE  CONFIGURATION  FOR  GENERATION  OF  AUTOCORRELATION 

SEQUENCE 


The  generation  of  a real  pulse-burst  code  which  will  have  a desired 
autocorrelation  behavior  has  then  been  equated  to  that  of  assigning  values 
to  the  coefficients  of  a nonrecursive  operator  so  that  its  unit-impulse 
response  best  matches  a desired  response  (i.e.,  the  ideal  autocorrelation 
sequence  (4.2-1)).  This  equivalency  is  significant  since  it  is  now  possible 
to  utilize  powerful  systems  theory  modeling  algorithms  for  generating 
excellent  real  pulse-burst  codes.  A particularly  effective  identification 
algorithm  for  this  purpose  will  be  presented  in  the  next  section. 

4.4  LINEARIZATION  ALGORITHM 

In  this  section,  an  algorithmic  procedure  for  selecting  the  coeffi- 
cients which  characterize  the  nonrecursive  operator  depicted  in 

Figure  4.3-1  so  that  its  unit-impulse  response  best  approximates  the  ideal 
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autocorrelation  sequence  R^n)  * E<$(n)  will  be  presented.  The  criterion 
to  be  minimized  is  squared  error  functional  (4.2-3)  in  which  the  element 
Pq  is  held  constant  at  some  prespecified  value.  The  ideal  code's  energy 
will  be  normalized  at  E = L+l  with  the  objective  of  encouraging  each  of 
these  conjectured  minimizing  code  elements  to  have  magnitudes  close  to  one 
for  reasons  of  efficient  signal  energy  transmission. 

The  method  to  be  used  is  the  so-called  linearization  algorithm  which 
has  been  found  to  be  particularly  effective  in  system  identification  and 
modeling  applications  [3]  - [4],  In  the  linearization  algorithm,  the 
signals  to  be  modeled  are  considered  as  approriately  defined  vectors. 

With  this  in  mind,  let  us  now  define  two  L+l  dimensional  vectors 


y[y]  = [y(0) ,y(l) ,y(2) , . . . ,y(L) ] 
y±  = [L+l,  0,  0,  ...  , 0]' 


(4.4-1) 


where  the  prime  denotes  vector  transposition.  In  these  representations, 
the  vector  y[p]  has  components  y(k)  which  correspond  to  the  causal  portion 
of  the  unit-impulse  response  of  the  given  nonrecursive  operator  and  there- 
fore the  causal  portion  of  the  autocorrelation  sequence  associated  with 
the  real  pulse-burst  code  {p^K  We  have  here  suggestively  expressed  the 
response  vector  as  an  explicit  function  of  the  Lxl  parameter  vector 

/ 

U ■ [w1.y2»  • • • (4.4-2) 

to  emphasize  the  dependency  of  y[p]  on  p.  The  vector  y^^  is  seen  to  be 
identical  to  the  causal  portion  of  the  normalized  ideal  autocorrelation 
sequence  that  is  to  be  approximated.  It  is  only  necessary  to  consider 
the  causal  portions  of  each  sequence  since  the  autocorrelation  sequences 
corresponding  to  real  pulse-burst  codes  are  known  to  be  even. 

The  functional  to  be  minimized  (4.2-2)  can  now  be  compactly  expressed 
in  the  standard- inner  product  format 


f(u)  - (y lul-yt>  W(y lwl-y1) 


(4.4-3) 


where  W is  the  (L+l)x(L+l)  diagonal  weighting  matrix  whose  diagonal  elements 
are  given  by  w(0) , 2w(l),  2w(2),  ....  2w(L) . It  is  noted  that  y l P ] — y 4 
is  the  error  vector  which  measures  the  amount  by  which  the  actual  code's 
autocorrelation  differs  from  the  ideal  autocorrelation. 

In  any  descent  algorithm,  the  existing  parameter  vector  is  perturbed 
in  a direction  so  as  to  cause  the  functional  being  minimized  to  take  on  a 
smaller  value.  With  this  in  mind,  let  us  consider  the  perturbed  functional 
value 


f(p+A)  = (y[y+A]  - y^  W(y[p+A]  - y^  (4.4-4) 

where  A is  a Lxl  parameter  perturbation  vector.  The  linearization  al- 
gorithm is  based  on  making  the  following  linearization  of  the  response 
vector  y(y+A) 

L 

y[p+A]  - y[p]  + l A . ay[p]  /3y  (4.4-5) 

k-1  * * 

where  the  A^  are  the  components  of  the  perturbation  vector  A and  the 
(L+l)xl  vectors  3 y[y]/i  are  the  so-called  "sensitivity  vectors”  of  the 
associated  nonrecursive  operator.  These  sensitivity  vectors  yield  the 
differential  amount  by  which  the  given  nonrecursive  operator's  response 
changes  when  a differential  variation  is  made  in  the  component  With 

reference  to  Figure  4.3-1,  it  can  be  shown  that  the  causal  portion  of  the 
sensitivity  vectors  are  given  by 

0y[vi]/5pk)  - [yfc,  pk_lt  ....  yQ,  0,  0,  ....  0]  (4.4-6) 

+ Iuk*  Pk+1 PL,  0,  0,  ....  0] 

for  k = 1,  2,  . . . , L.  The  determination  of  these  sensitivity  vectors 
for  a given  set  of  coefficients  (pq,  y^,  . . . , y^J  is  then  straight- 
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forward  requires  only  a simple  two  vector  addition  operation. 

If  the  lineariztion  of  the  perturbed  response  vector  as  given  by  ex- 
pression (4.4-5)  is  substituted  into  functional  (4.4-4),  a quadratic  (in 
A)  approximation  of  f(p+A)  results.  It  is  readily  shown  that  an  optimum 
selection  for  the  parameter  vector  perturbation  A to  minimize  this  quadratic 
approximation  must  satisfy  the  following  consistent  system  of  L linear 
equations  in  L unknowns 

[t'wt  1A  - T'wiy.-yM}  (4.4-7) 

u u yi 

where  T is  a (L+l)xL  matrix  whose  kC^  column  vector  corresponds  to  the 
sensitivity  vector  <*ylp]/3pk  for  k=l,  2,  ....  1.. 

It  has  been  shown  that  if  the  parameter  vector  p does  not  render  f(p) 
a relative  minimum,  then  the  pertirbed  vector  p+aA  will  yield  a smaller 
value  to  the  functional  being  minimized.  Here  a is  a step  size  scalar 
selected  in  the  interval  0 < a < 1 and  A is  a solution  to  relationship 
(4.4-7).  With  this  in  mind,  the  linearization  algorithm  as  described  in 
Table  4.4-1  is  offered.  This  algorithm  has  guaranteed  convergence  and 
in  fact  converges  in  a quadratic  manner  (e.g.,  see  ref.  [3]).  A computer 
program  to  implement  this  algorithm  is  given  in  section  8.2. 


_ " • • ' v.  . 
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Step  1 
Step  2 
Step  3 
Step  4 


Enter  initial  guess  y 
Calculate  9y[y]/l  y^  k = 1,2,.  . 

Evaluate  T 

y 

2 

If  ||  T7  W{y^  -y [ y ] } ||  is  less  than  some  prescribed  small  positive 

values  (e.g.,  10  ^) , then  the  algorithm  is  said  to  have  converged 

* 


otherwise  go  to  step  5 
Step  5 Solve  relationship  (4.4-7)  for  A 

Step  6 Sequentually  set  a = 1,  1/2,  1/4,  1/8,  . . . and  evaluate  f(y+aA) 

until  condition  f(y+aA)  < f(y)  is  first  satisfied 
Step  7 Let  y = y+aA  and  go  to  step  2 

TABLE  4.4-1  BASIC  STEPS  OF  LINEARIZATION  ALGORITHM 


4.5  INITIAL  CODE  SELECTION 

The  hypothesized  functional  f(y)  which  is  to  be  minimized  is  a quartic 
function  of  the  code  elements  y^.  One  would  therefore  anticipate  that  there 
will  exist  a number  of  different  code  selections  which  will  render  f(y)  a 
relative  minimum  (i.e.,  values  of  the  code  elements  U|,|ir),  . . . , y for 
which  the  gradient  of  f(y)  is  zero).  Some  of  these  "minimizing"  codes  will 
have  an  associated  autocorrelation  sequence  which  has  a desired  impulse- 
like character isitic  while  other  minimizing  codes  will  not  be  so  character  iz 
The  particular  relative  minimum  code  to  which  the  linearization  algorithm 
converges  is  to  a large  extent  dependent  on  the  initial  code  selection 


0 

1* 


} 


The  quantity  T W{y.-y[y]}  is  in  fact  equal  to  minus 
y i 2 

of  f (y ) (i.e. , -0.5V^f  (y) ) while  the  notation  ||  x || 


one  half  the  gradient 
stands  for  x/x. 
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which  is  made  at  step  1 of  the  algorithm  as  specified  in  Table  4.4-1. 

Since  the  linearization  algorithm  is  of  a descent  nature,  it  follows  that 
when  the  initial  code  selection  has  a satisfactory  autocorrelation  behavior, 
then  the  code  to  which  the  algorithm  converges  must  be  of  a superior  quality. 
Thus  the  initial  code  selection  is  the  most  critical  factor  in  generating 
suitably  good  pseudo-optimum  codes. 

For  code  lengths  in  which  there  exists  a Barker  code,  a natural  initial 
code  selection  would  be  the  Barker  code.  On  the  otherhand,  an  appropriately 
chosen  real  Huffman  code  would  suffice  for  any  code  length  (e.g. , see  ref. 
[1]).  Alternatively,  a concatenation  of  shorter  length  Barker  codes  has 
also  been  found  to  be  effective  in  constructing  reasonably  well-behaved 
initial  codes  of  long  length. 

In  formulating  the  functional  to  be  minimized,  the  code  element  p^ 
was  held  fixed  so  as  to  prevent  the  algorithm  from  converging  to  an  unde- 
sired optimum  code  of  length  one.  It  has  been  empirically  determined  that 
different  choices  of  p^  will  result  in  codes  which  has  significantly  dif- 
ferent energy  efficiency  factors,  normalized  maximum  side  lobe  ratios, 
and,  code  variances.  With  this  in  mind,  the  linearization  algorithm  is 
run  with  a number  of  different  choices  for  p^  taken  from  the  set 

{6,  26,  . . . , p6 } 

where  typically  p6=l  with  6 being  a positive  number  less  than  one  (e.g., 
5=0.1).  From  the  family  of  p pseudo-optimum  codes  which  results  from 
each  of  these  fixed  p selections  of  p^,  it  is  generally  possible  to 
select  one  which  is  deemed  superior  (e.g.,  one  whose  energy  efficiency 
is  a maximum) . 
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SECTION  S 


NUMERICAL  RESULTS 

5.1  PSEUDO-OPTIMUM  CODE  OF  LENCTH  11 

The  linearization  algorithm  as  described  in  Section  4.4  was  first 
used  to  generate  a pseudo-optimum  code  of  length  11.  It  is  noted  that  for 
this  code  length,  there  exists  an  associated  Barker  code.  With  reference 
to  Table  3.4-1,  a logical  choice  for  the  linearization  algorithm's  initial 
code  is  then  given  by 

P°  = [V°,  1,  1,  -1,  -1,  -1,  1,  -1,  -1,  1,  -1)  (5.1-1) 

where  it  is  noted  that  the  Barker  code's  initial  element  value  1 has  been 

o 

replaced  by  the  fixed  parameter  p^.  The  linearization  algorithm  was  then 
used  to  generate  the  pseudo-optimum  codes  which  corresponded  to  each  of  tin1 
following  ten  selections  for  p”. 

Pq  = 10.1,  0.2,  0.3 1} 

In  this  algorithmic  procedure,  the  error  weights  w(n)  were  all  set  to  one. 
The  normalized  maximum  side  lobe  ratio,  the  energy  efficiency  factor,  and, 

the  code  variance  (as  defined  in  Section  3.3)  which  corresponds  to  each  of 

o 

the  pseudo-optimum  codes  iteratively  generated  for  these  p^  selections  are 

shown  in  Table  5.1-1.  The  best  code  from  this  family  is  deemed  to  be  the 

one  associated  with  p”  = 0.5  since  maximum  energy  efficiency  utilization 

° 1 

was  of  primary  concern.  A plot  of  this  pseudo-optimum  code  and  its  as- 
sociated autocorrelation  function  is  shown  in  Figure  5.1-1. 


lr  maximum  side  lobe  reduction  were  of  primary  concern,  then  the  code 
corresponding  to  Mq  = 0.1  would  have  been  chosen. 
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TABLE  5.1-1  INITIAL  CODE  SELECTION-LENGTH  11 


In  terms  of  normalized  maximum  side  lobe  ratio,  the  pseudo-optimum  code 
shown  in  Figure  5.1-1  is  superior  by  6.36db  over  the  Barker  code  of  length 
11,  and,  by  6.00d.b  over  the  complex  Huffman  code  of  length  11  as  shown  in 
Figure  3.5-2.  Although  this  is  but  one  example  of  a short  length  real 
pulse-burst  code,  it  does  illustrate  how  one  may  generate  effective  real  codes 
using  an  algorithmic  approach. 

To  further  compare  the  inherent  signal  characteristics  of  the  pseudo- 
optimum, Barker,  and  Huffman  codes  of  length  11,  their  associated  ambiguity 
functions  were  calculated  by  means  of  relationship  (3.1-5)  and  are  shown 
in  Figure  5.1-2.  In  these  plots,  the  ambiguity  surface  is  sliced  along  the 
w=0  axis  so  that  the  associated  autocorrelation  function  corresponds  to  the 
leading  edge  of  the  surface  as  is  readily  apparent  in  the  cases  of  the  Barker 
and  Huffman  codes.  A careful  examination  of  these  ambiguity  surfaces  reveals 
that  the  pseudo'-opt imum  and  Rarker  codes  are  clearly  superior  to  the  llulTm.m 
code.  More  over,  the  pseudo-optimum  is  somewhat  better  behaved  than  the 
Barker,  particularly  relative  to  autocorrelation  behavior. 
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5.2  PSKUDO-OPT  1 MUM  CODE  OK  LENGTH  128 

The  true  effectiveness  of  the  propi>seti  algorithm  for  synthesizing 
real  pulse-burst  codes  Is  made  evident  when  constructing  long  length 
codes.  To  demonstrate  this  assertion,  a real  pulse-burst  code  of  length 
128  was  next  generated.  When  considering  this  relatively  long  length 
code,  the  choice  of  the  initial  code  p^,  p^,  . . . , y®  used  in 

step  1 of  the  algorithm  becomes  even  more  critical  as  well  as  difficult 
to  make.  In  this  example,  the  real  code  of  length  128  as  tabulated  in 
reference  II]  was  used  to  initiate  the  algorithm.  This  initial  code 
has  excellent  characteristics  as  shown  in  Table  5.2-1. 


Normalized  Maximum 
Side  Lobe 

Energy  Efficiency 

Code  Variance 

X 

P 

F 

Initial  Code 

9.508  x 10~4 

0.3348 

0.2757 

JLseudo-Opt  imum 
Code 

— ft 

8.125  x 10 

0.3426 

0.2862 

TABLE  5.2-1  RELEVANT  CHARACTERISTICS  OF  REAL  PULSE-BURST  CODES  OF  LENGTH 

128 


Using  this  initial  code,  the  algorithm  was  found  to  converge  to  a 

pseudo-optimum  code  whose  normalized  maximum  side  lobe  ratio  was  found 

-4 

to  have  decreased  by  a factor  of  10  over  that  of  the  initial  code. 

The  relevant  parameters  which  characterize  this  pseudo-optimum  code  are 
also  shown  in  Table  5.2-1.  In  addition,  a plot  of  the  codes  magnitude 
and  its  associated  autocorrelation  and  ambiguity  functions  are  given  in 
Figure  5.2-1.  It  must  be  mentioned  that  these  excellent  characteristics 
and  ambiguity  function  behavior  can  only  be  approximated  in  practice 
due  to  the  inherent  incapability  of  exactly  reproducing  the  precise 
amplitudes  for  the  code  length  required. 
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(a) 


(b) 

Ruw 
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! 

-t 
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-128T 


0 


128T 


Figure  5.2-1  Pseudo-Optimum  Code  of  Length  128 

(a)  Magnitude  of  Code  |p(t)|,  (b)  Corresponding 
Autocorrelation  Function  R (t) 
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SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  iterative  procedure  as  herein  described  for  synthesizing  real 
pulse-burst  codes  has  proven  to  be  effective  in  all  examples  treated. 
Namely,  the  pseudo-optimum  code  to  which  the  algorithm  converges  has  been 
found  to  posses  such  desirable  characteristics  as  having  an  associated 
impulse-like  autocorrelation  function  and  a well-behaved  ambiguity  func- 
tion. Furthermore,  the  algorithm's  convergence  rate  is  quadratic  in 
nature  thereby  offering  a computationally  efficient  code  synthesis 
procedure.  This  procedure  therefore  promises  to  be  a very  useful  too] 
for  constructing  good  real  pulse-burst  signals. 

To  achieve  the  algorithm's  true  potential,  however,  it  is  essential 
that  a systematic  method  for  constructing  a suitable  initial  code  be  de- 
veloped. The  method  suggested  by  Akroyd  [1]  appears  to  satisfy  this  re- 
quirement although  further  investigations  along  this  line  need  be  made. 

For  example,  it  might  be  possible  to  further  adapt  some  of  the  techniques 
as  given  by  Schroeder  [7]  to  construct  appropriate  initial  codes.  The 
importance  of  this  initial  code  selection  cannot  be  over  emphasized. 

Another  consideration  which  must  be  made  is  that  of  generating  codes 
whose  energy  efficiency  is  appropriately  large.  Initial  investigations 
of  how  one  might  achieve  this  capability,  while  still  maintaining  the 
algorithm's  basic  features,  have  been  made.  In  essence,  this  requires 
that  an  inner  product  which  gives  a positive  measure  of  energy  efficiency 
be  added  to  the  functional  f(p).  A minimization  of  this  new  functional 
would  then  tend  to  simultaneously  generate  a real  pulse-burst  code  which 
has  a relatively  large  energy  efficiency  utilization  as  well  as  an  impulse- 
like autocorrelation  behavior. 
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SECTION  S 


8.1  APPENDIX  A 

AMBIGUITY  FUNCTION  PLOT  PROGRAM 


1 


C 

c 


c 

c 

c 

v- 

c 

c 


PROGRAM  MAIN(INPUT» OUTPUT  *TAPE5-InPJT 
i , r APE6-0JH JT 

t»TAPEl*TAPE2»TAPE3»TAPE4»TAPE7»TAPEd*TAPt9) 

**♦  THIS  PROGRAM  GENcRATES  He  AMBIGUITY  ♦ 
*♦*  FUNCTION  AND  PLOTS  IN  THREE  * 

***  DIMENSIONAL  SPACE.  * 


COMPLEX  0»dlC0M,B2C0M»0rtMP.XTcMP,E<l,EX2#CSUM 

i > A » T t MP 

DIMENSION  D(i.23),dlC0M(12  5)»d2CJMC123J#X(2  5t»l) 
i »MASK(A001),0aTA(236) 

* »0IN(  128>,TtMP(  43),A(43»128) 


REWIND  1 
REWIND  2 
REWIND  3 
REwIND  A 
RtwIND  7 
REWIND  8 
REWIND  9 

PI-3. 141 592oS3SdR8 


♦**  T 

L 

♦♦♦  NF 

♦ * ♦ 

* * ♦ 

***  DIN  — 
♦ ♦ * 

♦ 0 


SAMPLING  TIME 
NUMBER  OF  CODE  LENGTH 
NUMBER  OF  LINeS  FOR  EVtRY 
2 *P I / T WHICH  ARE  PLUTTsD 
PARALLEL  TO  HE  TIMl  AXIS 

real  part  of  each  code 

t L E M E N T 
CODE  ARRAY 


T-l. 


L»1  23 
AL  -L 

LI-L/3+1 
L ? ■ L 1 * 2 
L h»L  / 2 

L 2 M» ( L-l ) *2-1 


NF  *fc  A 
NF2«NF  *2 
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N F 3 • N (■  ♦ \ 
KF  c^NF 
NF  Mi »NF-1 


K|  A0(  2» ♦ I (OlN(i)#M»l»L) 

0 J 14  M ■ l * l 

14  0<M)*CMPLX(QIN(M)  ,0.0) 

CSUM-CMPLX ( 0.0, 0.0) 

DT  30  M«l,L 

ij  ccu**csu,ucnNjG<o(i)  >*D(  d 
SU»«RE AL (CSUrt) 

SC4-SORHAL/SUM) 

DQ  41  M ■ 1 » L 
4j.  D<M)»SCA*D(1) 

• RITE (6, LOO) 

loo  FORMAT (/, 3 X, ♦ CODE  0 *) 

^IM  * , (0(«),M»1,L) 

O£»0. 

***  S u ^ A T I IN  INVOLVED  IN  A(i3  IGJl  T Y 

***  function  is  Computed  a-'  fft. 

00  7 N*l>2 
0 0 1 I P « 1 , L 
CALL  StCONO(TIMt) 

PJIM*,"TI^E  IS  ",Tm 
1 1 * 1 p- 1 

NNF  *2* ( I P+ ONc ) 

IF (NNF.LO.o)  oJ  TO  21 
00  2 M ■ 1 , N N F 

2 OATA<M)*CN?LX <0. 0,0.0) 

?i  I S * I P ♦ONE 

IF(IS.GT.L)  GO  TO  22 
00  3 N • I S » L 
*T*tf*2-H-2 

DTcMP*0(N)+CoNJG(0( 1-II-jNe ) ) 

■j  DATA(MT)«R£AL(DTtMP> 
DATA(MT*1)*AIvIAG(DTEMP) 

3 CONTINUE 

22  LL«2*»l+l)*l 
DO  Id  M«LL,NF£ 

DATA(M)*0.0 
Id  CONTINUE 

IF (L.LE.NF  ) GO  T J 2 3 

COMPUTES  8Y  NOT  JSlNS  FFT. 

27  CALL  AFOUP<OATA#NF,L#uNt, T) 

GO  TO  2b 

l)  CALL  FOUR  i ( 0 fi  T A,  NF  » -1 ) 

2 o IF ( IP.GT.L1)  GJ  TO  31 
DL  6 N»l,NF 
"T*M2-1 

6 A(IP,M>«CMPLX(DATA(1T).UATA(MT«-1>) 

IF (IP.Nfc.Ll ) GO  TO  1 
DO  3?  MJ ■ 1 * NF 
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■ 


DJ  33  MI-1, LI 
33  TfcNP(MI)«A<MI,MJ) 

IF(K.EQ.2)  GO  TC  3b 
WRITE! 1 * + ) TEMP 

GO  T 0 32 

3b  wRITE(9,*)  TcMP 
32  CONTINUE 
GD  TO  1 

31  IF(IP.GT.L2)  GO  TO  <,3 
1 1 P » I P-Ll 
CO  9 M«1 ,NF 
MT«F*2-1 

9 A(IIP,M)»CMPLX(OATA(MTi,OATAlMT+i)> 
IF(IIP.NE.Ll)  GJ  TO  1 
DO  39  MJ«1,NF 
DO  35  Ml  * 1 » L 1 
35  TEMP(MI)«A(MI,MJ) 

IF (K.EQ.2 ) GO  TO  37 
WRITt(3,*>  TEMP 
GU  TO  39 

37  w»ITE(7,+)  TEMP 
39  CONTINUE 
GO  TO  1 
9 3 I IP- IP-12 
CO  99  M«1,NF 
MT-M+2-1 

99  A(IIP,M)-CMPlX(DATA(MT),DATA(MT*1>) 

IF(IP.NE.L)  GO  TO  i 
DO  95  M J - 1 , N F 
DO  96  MI-1, LI 
96  TE*P(MI  )-A(MI,MJ) 

IF  ( K . E 0.  2 > GU  TO  99 
WR I TE ( 8 , * ) TEMP 
GO  TO  95 

99  FR  I T t ( 9,  * ) TEMP 
96  CONTINUE 

1 CONTINUE 
ONE  «i . 

7 CONTINUE 

C *+*  THIS  PROGRAM  HAS  THE  CAPABILITY  JF  * 

C PLOTTING  FOR  T Ht  INTERVAL  -6*PI/T  TU  ♦ 

C *+*  ♦b+PI/T  ALONG  OMEGA  AXIS.  ♦ 

CALL  PLOTS 

I STEP-10 

ISTFPl-ISTEP+l 

DEL-T/ISTcP 

IPOINT«L*ISTlP 

L2M»( IP0INT*2+1  )/10 

IP21-IPOINT»2+i 

I6AS6-NF3 

NF3P-NF3+1 

NF6P-NF6+1 

NFD-NF+5+1 

IL INE-0 

CALL  PLOT3D(10,O.J,X,3.u,G.JJ23,u.w3125 
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f > o 


i #-u.  J'3o»  IllNtf  IP21 

f #-<*3. 0» -2  j.  J»  10.0*  . 25  >20.0,. MASKED) 

P t w I NO  i 
REwINO  3 
REMIND  A 
REWIND  7 
REWIND  8 
REWIND  9 

C EVERY  LINE  IS  DKAwN  PARALLEL  TO  Trlt  * 

***  TINE  AXIS.  THE  LINES  WHlCri  SHJULD  3c  * 
♦**  HI  ODEN  BY  PREVIOUS  LINES  ARE  MASKED.  ♦ 

DO  13  I-NF3P.NFD 
IL INE-I-NF3 
TR "NFfcP+I-I 
II«NF3*1-I 
1 1 «-I I 

0MG»PI*II*2./(NF*T) 

EXl«CMPLX(J.j,0MG*T/2. ) 

Fx2»CMPLX<0.0,0MG*T  ) 

IF(IR.LE.NF3)  GO  TO  16 
IND-IR-IBASE 
IF(IND.NE.i)  GO  TO  12 
I F ( I , EQ. NF3P ) GO  TO  12 
IBASE-IBASE-NF 
GU  TO  12 
lo  IND»I-I3ASt 

IF (INO.NE.NF ) GO  TO  12 
IBASE*IBASE+nF 
12  IF { INO.NE.l ) GO  TO  50 
PEwIND  1 
REWIND  3 
R t W I ND  A 
REWIND  7 
REWIND  9 
REWIND  9 

-j')  R E AD  ( 1#  ♦ ) TEMP 
DO  33  M ■ 1 » L 1 
33  31C0M(M)«TEMP(M) 

HEAD(3t*)  TEMP 
DO  39  M* 1 » L 1 
MT«L  UM 

39  31C0M(MT)«TEMP(M) 

R E AD ( 0 » + ) TEMP 
DO  A 7 M«1»L 1 
MT-L2+M 

A 7 B1C0M(MT)«TEMP(M) 

READ ( A.*  ) TEMP 
DO  A 2 M*1,L1 
A2  B2C0M(M)«TEMP(M) 

Rt AD ( 7»  + ) TEMP 
DO  AO  M»1,L1 
MT»L1+M 

AD  62CCM(nT)*TEMP(M) 

READ ( 9» ♦ ) TEMP 
DO  A3  M • 1 » L I 

50 





M » L 2*M 

43  82C0m(MT)»TEMP(M) 

xTLMP«T  + C£XP(EXl)*SINC(G'1G*T/2.>*diCJMm 
X<IPDINT*1)»CA9S(XTEMP> 

DO  10  J«1»L 
J J » J- 1 

DO  10  <»1»  ISTEP 
*K«K 

EXl»CNPL  X (O.C»OMG+(  T-DEL*i<i<)  / 2#  ) 

E X 2»C  MPL  X(0.0,0MG*<  2.  ♦T-OcL*f<K)  ✓ 2 • > 
XTEMP-(T-DEL*KK)*CEXP(cXi) 

1 +3INC(OMG*(T-D£L*KK)/2»)*31CGM(J) 

* ♦ DEL*KK*CEXP(EX2) 

t *SINC  (0MG*0EL*«K/2.  )#B2C0M(  J ) 
*X*CABSIXTEMP) 

KT« J J+ISTEP+K 
X( I POINT+l-KT ) «XX 
X( IP0INT+1*KT)«XX 
10  CONTINUE 

CALL  PL0T3D! 10»C.O»X,0.0,G.Gu23,0.03125 

* »-0.036,ILINE»IP21 

J »-<t5.0,-20.C»10.0».25,2w.J,MASK,C) 

IF (I .NE.MF3P)  GO  TO  13 
24  WRITE (6» 140) 

140  FORMAT!/, 3X, ♦OMEGA  » J* 1 
DO  23  IM»1,L2M 
BASfc»10+( IM-1 ) 

IBS»3ASE*1 

I8L»3ASE*10 

PRINT  ♦, (X!M),M«I3S»I8L) 

23  CONTINUE 
13  CONTINUE 
CALL  cFPLOT 
STOP 
END 


SU3R0UT  IM€  AF0US{DaTA»NF»l#3NE»  T) 


***  COMPUTES  FQjSIcR  TRANSFORM  * 

DIMENSION  0ATA(256)»DT ( 2 56  ) 

COMPLEX  t X » S J M * T M P 
P I » 3 . 1A15326535898 
00  1 I P» 1 * NF 
I P P»I p *2-1 
1 1 P - 1 P-1 

SUM»CMPLX(0. 0*0.0) 

L P «L  * 1 
DO  2 N»2,LP 
1 1 « 2*N— 1 
NN-K’-l 

X«-2.*PI*IIP+NN*T/NF 
TMP*CMPLX (0.0*  X ) 
tX-CEXP(TMP) 

TMP-CMPLX(0ATA(II)»DATA(II*1)  ) 

SUM*TMP*fc  X*SUM 

2 CONTINUE 

DT ( IPP )»«cAL(SUM) 

DT(IPP*1)»AIMAG(SU'1) 

1 CONTINUE 
NF2  *NF  *2 
DO  3 I * 1 » NF  2 

3 D ATA ( I ) »DT  < I ) 

TURN 

END 
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8.2  APPENDIX  B 


PSEUDO-OPTIMUM  CODE  ALGORITHM  PROGRAM 


PROGRAM  MAIN(INPUT»0UTPUT»TAPc3*lNPUT 
i , lAPtfc-OJTPUT 

4 » TAPei. TAPe2» TAPE3 ) 

***  THIS  PROGkAM  CGMPUTcS  T HE  OPTIMUM 
♦♦♦  PULSE  BURST  CODE. 

COMMON  R(255),N.NN»AN»0(235)»0J»NP 
01  ME  NS  I ON  D(i27).  00  (i27)»P(i27),0T(i.27) 

5 »S(127),30(255).DIN(127) 

REAL  LAM 


+ **  D — 
***  R — 
*♦*  N — 
* * t 

***  0 -- 
*** 

♦**  LAM 


♦++  BET 
**♦ 


CODE  ARRAY 

CORRELATION  ERRROR  ARRAY 
THE  LARGEST  INDEX  OF  CODE 
THE  LENGTH  OF  CODE  IS  N*1 
WEIGHTING  FACTOR  OF  MEAN 
SQUAR  ERROR 

ADJUSTABLE  FACTOR  TO  GET 
BETTER  EFFICIENCY 
(WEIGHTING  OF  PENALTY 
FUNCTIONAL ) 

ADJUSTABLE  FACTOR  TO  GET 
BETTER  SPKt  TO  SIDtRJBc 
RATIO 


N-127 
NP«N*1 
nP2«N«-2 
NH»N/2 
NHK»NH— 1 
AN  *N 

NN«  2*N  + 1 
L A M » 0 • 0 
3FT«50. 


WR I T£ ( b»  35C ) LAM 

350  FO»MAT( /»3X,*LAMBDA  « *#Fi0.5) 
WRITt(6»360)  BET 

360  FORMAT(/»3X#+8ETA  » *»Fi0.5) 
3*»U  FORMAT  ( / > 3X  » ♦ J MATRIX  s +) 
READ*, (3(M)»M*1#2I) 


« E AD ♦ # (3(1)  ,W2»h  1) 

03  32  M«R2,nn 
32  0 ( M ) * 1 • 

00  35  M-l.NN 
35  CG ( M ) ■ Q ( M ) 

READ  *»00» 101M) * M » 1 » M » 

SUM«  OC  *DG 
00  60  M«1,N 
6C  SUM»SUM*0<M)*0< M> 

F A C « ( AN  + 1. )/SUM 
F AC*SORT<  FAC  ) 

DO  *00+  FAC 
00  61  M«1»N 
6 i D(M)«D<M)*FAC 

PER0»C0RRE(0,S*LAM) 

CALL  PVEC(P»GRA0»0»LAM*S) 

1-0 

CALL  CONVT ( R»  N ) 

WRITE<6,100>  I»PERG,GRAD 
WRITc (6*200) 

CALL  WR(P>M 
WRITE ( 6 » 300) 

PRIM  *,00 
CALL  WR(0,N) 

WP ITE l 6»  310 ) 

CALL  WR ( R » NN ) 

CALL  EFVAOQ,D»R(M*l>,N) 

CALL  CQNVT(«,M) 
on  20  M*1»N 
20  OD (M>-D«M> 

I-l 

Ml»l 

♦**  COMPUTES  P VECTOR 

ll  CALL  PVF.  C(P»GRAD*D»LAM»S) 
PER0»C3RRt<D. S»  LAM) 

00  66  M«1>N 
66  SUM»SUMP(M)*P(M) 

ALFMAN+1. ) /SUM 
ALF*SQRT( ALF  ) 

IF (ALF.GT.2. ) A LF  » 2 • 0 

SEARCHING  THE  STEP  SIZE 

68  00  13  J»1»1C 
PRINT  *>ALF 
DO  1A  K * 1»  N 

IR  0<K)»DD(K>*ALF*P(K> 
ER0»C0RRE(D,S»LAM) 
IF(ERO.LE.PfcRO)  GO  TO  15 
ALF-ALF/2. 

13  CONTINUE 
15  PERu-ERO 
DO  18  M » 1 » N 
13  D0(*)»0(M) 

SUM«0. 


A L L « ( 1./2. ) **R 
IF  (ALF.LE.ALL  ) GJ  TJ  17 
IF(GRAO.LT..OOC1)  SO  TO  17 
IF ( J. tO. 10)  GO  TO  17 
IF (I.EG.50)  GJ  TO  17 
IM*1 
GO  TO  11 

♦**  AO J OS  T S WEIGHTING  J 

17  SUM*0. 

00  33  M«1,NN 

33  $UM«$UM*R(M)*R(M) 

00  3 ^ M»1,NN 

34  Q(M)*QQ(f()*6£T*R(M)*R(i't)/SUM 
*R I TE ( fe» 100 ) I » PE  R0»  GR AG 

100  FORMAT  (//,10X»*ITER.  - *,I3,5X, 
i * ERROR  » *.£11.3, 

3t>X  » *GR  ADI  E NT  » *,F20.S) 

WRITE (6.20C ) 

20C;  FjrmaT  (IX.  * P VECTOR  :♦> 

CALL  w R ( P » N ) 

WRITE (6,300) 

300  FORMAT!/, IX,*  CODE  0 :♦) 

PRINT  *,00 
CAlL  WR(O.N) 

CALL  C ON V T l R , N ) 

WRITE(6,310) 

310  FORMAT!/,  Ut*  CORRELATION  R 
CALL  WR(R»NN) 

CALL  EFVA (D0»0»R(N*1 >,N) 

(•RITE  (6,340) 

CALL  «R(3»NN> 

♦**  COMPUTES  CORRELATION  FOR  TRUNCATED 

**♦  CODE 

CALL  T R UNC ( DO » 0 , N ) 

P T R A« CQRRE ( D, S»  LAM) 

CALL  CJNVT(R.N) 

«RITE(b,3lO) 

CALL  wR(R,NN) 

CALL  CONVT(R.N) 

IF(K1.Ej.4)  GO  TO  38 
IFIM1.E3.1)  GO  TO  38 
Ml ■ M l ♦ 1 
GO  TO  11 

3a  STOP 
ENO 


* 

* 


Ct  o 


SU&RUUTINE  PVeCtP.GRAD.D.LAM.S) 


♦ + ♦ COMPJTES  F R EC  He.  T MAWIX  ♦ 

*♦*  AND  FINOS  P VECTOR.  * 

COMMON  R(255)»N#NN»AN.J(255)»J0»NP 

01  MENS  ION  FRM(1  27.  127)»P(127).$U27), 3(127) 

$ » D ( 1 27 ) 

PEAL  LAM 

DO  1 J ■ 1 » N 

US»N*1-J 

U L ■ U S ♦ N 

DS*  J*1 

DL  * D S*N 

DO  2 I ■ 1 » NN 

IF  (I.GT.UL)  GO  TO  3 

III- I -US 

IF(III)  3, A, 5 

3 UK*e. 

GO  TO  6 
A UK-DO 
GO  TO  6 
5 1 1 - I -US 
UK»D( II ) 

b IF(I.LT.DS)  GO  TO  9 
III-I-DL 
IF(III)  7.8.9 
7 1 1 — 0 L — I 
OK -0  < II) 

GO  TO  10 
3 DK  «D0 
GO  TO  10 
9 OK  «0 • 

1C  IF(I-NP)  30.39.2 
33  F R M ( I . J ) * UK ♦DK 
GO  TO  2 
39  S(J)«UK*DK 
2 CONTINUE 
1 CONTINUE 
REWIND  1 
DO  23  I«1.N 
DO  2 A J • 1 # N 
SUM«Q. 

DO  25  K » 1 . N 

25  SUM»1,*FRM(K»I)AFRM(K.J)*J(K)+SUM 
PIJ)«Sl)H*8(N*l)*S(I)*S(J) 

2A  CONTINUE 
wRITt  (1)  P 
23  CONTINUE 
DO  26  I ■ 1 » N 
S I J M « 0 . 

CO  27  K * 1 » N 

27  SUM«FRM(K»I)*R(K)*0(K)*1.+SUM 
B(I)»SUM+0(N+1)*R(N+1)*S(I) 

26  CONTINUE 
REWIND  1 

DO  36  I-l.N 

56 


;v.£.l.-g  


A 


READ  (1)  ? 

Du  3 7 J ■ 1 » N 
37  FRM  ( I , J)-P(J) 

3 1  CONTINUE 
OH  31  I » 1 # N 

31  FRM( I, I )-LAM  + FRM(  I*  I > 

SUM-O. 

OC  32  I « 1 » N 

32  SUM-SUM+D( I >*D(  I) 

SUM-SORT (SUM/AN) 

DO  33  I * 1 » N 

33  S(I)«SI6N(SUM,0(I); 

♦**  COMPUTES  GRADIENT 
DO  3A  I ■ 1 » N 

3 ^ 3(n*B(I)-LAM*(D(I)-S(in 
SUM-O. 

DO  13  I ■ i » N 
13  S U M » S UM ♦ 8 ( I ) * 3 ( I ) 

SUM-SUM*^. 

GR AD-SQRT(SJM) / AN 
CALL  S1MQ(FRM»B»N»  ING) 

I F ( ING.cO.O)  GO  TO  35 
WRITE  (f>»lJl) 

lul  F OR  MT(/»3X»*  MATRIX  IS  SINGULAR*) 
35  DO  17  I-i,N 
17  P(I)«8(I) 

RETURN 

END 
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<~>  o 


FUNCTION  CORRE <D.S>LAM) 


* + * COMPUTES  THE  CORRELATION  bF  CODE 
♦♦♦  AND  MEAN  SOUAR  ERROR. 

COMMON  RI253)»N>NN,AN»J12;>5)>00»NP 
DIMENSION  D(i27),S(127) 

REAL  LAM 
Nl-N-1 
TO  1 I»I#Nl 
1 1 «N+1*I 

R(II)»D(I)*Ou 

JJ«N*1-I-1 
DO  2 J»1»JJ 

2 Bill )>R ( II > *0( J )*DI J*I ) 

1 CONTINUE 

R IN*1)»DC+Dc 
0 0 3 I « 1 * N 

3 S ( N*1  ) «R<  N*1  H-D  ( I ) *D  ( I ) 

R ( NN ) * 0 ( N ) *00 

00  A I ■ 1 » N 
1 1 * NN* 1— I 
A R ( I) *R  (II  ) 

00  5 I * 1 » N 
1 1 *N+ 1 ♦ I 
R(I)«-R(I ) 

5 P(II)«-R(II) 

P (N  + l ) » AN-R (N*l >+l. 

F R U*0 . 

DO  b I • 1 » NN 

6 E R 0*  0 ( I ) *R  ( I ) *R  ( I ) *ERU 
S U M »0 • 

DO  7 I * 1 / N 

7 sum- sum* iDc i j-s (i n * ( o ( i i-sm ) 

ERO»ERO*LAM*SUM 

CQPRE-ERO 

RETURN 

END 
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bUBRJUTINE  cFVA(DO, D>RX# N) 


C ♦♦♦  COMPUTES  EFFICIENCY  AMD  CJDc  * 

C **♦  VARIANCE.  * 

DIMENSION  0(1) 

AN  *N 

D M ■ DO ♦ DO 
DO  2b  M«1,N 
OQ»C (M)*D(M» 

2b  PM»AMAX1(DM,D0) 

E F ■ R X / ( ( AN  + l , ) * DM  ) 

WR1TE16.320)  EF 

320  F0RMAT(/#1CX#*EFFICIENCY  » *»FiC.7) 

FF«BF*DM 

EF-SQRT(EF) 

SUM»(DO-EF ) **Z 
DO  26  M«1,n 
A90*A8S(D(M)  ) 

26  SUM*  ( A3D-EF  ) **2  + 3U.M 
VAR-SUM/RX 
WR I T E ( 6»  330  ) VAR 

330  FORMAT ( / »lUX»*C0D6  VARIANCE  * +>F10.7) 
kE  TORN 
END 


bUBRQUTINE  TR'JNC(Xi.»X»N) 

C **♦  TRUNCATcS  THE  CODE. 

DIMENSION  x (1 ) 

X1«X1*1C0. 

IX1«X1 

X1-IX1/100. 

DO  1 I ■ 1 * N 
X<I*-X(I)  ♦100. 

I X » X ( I ) 

i xm-ix/ioc. 

RETURN 

END 
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SUM  OUT  iNt  C UN V I ( R » N ) 


++*  CONVERT  TO  CORRELATION  FROM 
***  CORRELATION  ERROR  OR  VICt  VtRSA. 

DIMENSION  R ( 1 ) 

A N * N 

DO  23  M»1,N 
1 1 “2+N+2-M 
R ( M ) *-R ( M ) 

23  k ( II )»R(M) 

R ( N ♦ 1 >»AN-R < N + l )*1  . 

RETURN 

END 


SUbROUTlNt  WR ( X X»  NNN  ) 

**♦  THIS  SUbROUTlNt  IS  USED  TO  PRINT 
♦»*  OUT  THE  RESULTS  FOR  CONVENIENCE. 

DIMENSION  X X { 1 ) 

NL AS»NNN/10 
DC  bl  I«1,NIAS 

is«ic*<i-im 

I L * I S ♦ 9 

WR I TE ( 6»  3 70 ) ( X X ( J)» J»IS»  IL) 

37U  FORMAT ( / » 3X * 10 ( 1 1 1 . 2 ) ) 

5 i CONTINUE 
IIS-NLAS*10 

IF(IIS.EO.NNN)  GO  TO  I 
IS«NL AS*10+1 

WRITE(fe»370)  (XX(J), J»IS»NNN) 

1 RETURN 
END 
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